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A new simple model exhibiting a noise-induced ordering 
transition (NIOT) and a noise-induced disordering transition 
(NIDT), in which the noise is purely multiplicative, is pre- 
sented. Both transitions are found in two as well as in one 
dimension (where they had not been previously reported) . We 
show analytically and numerically that the critical behavior 
of these two transitions is described by the so called multi- 
plicative noise(MN) universality class. A computation of the 
set of critical exponents is presented in both d = 1, and d — 2 
(where they have not been previously measured). 

PACS: 05.40.+,] 



I. INTRODUCTION 

It is well established that the effects of noise in stochas- 
tic equations can be rather surprising and counterintu- 
itive. The presence of a noisy term can modify in a quite 
an unexpected way the behavior of the corresponding de- 
terministic (noiseless) equation. A well known example 
of this is the stochastic resonance in which the response 
to an external periodically oscillating field of a stochas- 
tic system is enhanced by increasing its internal noise 
amplitude 111. Other examples are the resonant acti- 
vation 0, noise-induced spatial patterns [pi, the noise 
induced ordering transition (NIOT) |4|-||l, and the noise 
induced disordering transition (NIDT) 0. This two last 
transitions were first analyzed together in an interesting 
work by Van den Broeck, Parrondo and Toral Q (VPT 
hereafter). They considered a Langevin equation with 
a monostable deterministic term, to which a non-trivial 
noise is added. Owing to the combined effect of the noise 
and the spatial coupling two real phase transitions were 
shown to appear: for small noise amplitudes the system is 
disordered, while by increasing the noise amplitude the 
system gets ordered in a symmetry-breaking state (ex- 
hibiting therefore a NIOT) . If the noise is increased fur- 
ther the ordered state is destroyed, (and therefore there 
is a NIDT). Afterwards, different works on NIOT and 
NIDT appearing in different models and systems have 
been proposed in the literature |Q-^. In particular, in [Q, 
a noise-induced reentrant transition is reported for a time 
dependent Ginzburg-Landau model with both additive 
and multiplicative noise. In MM first order transitions 



are shown to be induced by the presence of non-trivial 
noises. 

An interesting task is that of understanding the na- 
ture of this striking phenomena in the simplest possible 
model: understanding whether it comes from an inter- 
play between additive and multiplicative noises or not, 
and determining whether the critical behavior of these 
transitions can be described by any standard universality 
class. In this direction, the following Langevin equation 
has been proposed by VPT as a possible simplest model 
exhibiting these two noise induced transitions: 



d(j){x,t) 



dt 



= F[<j>{x, t)] + DW^(f> + G[(l){x, t)]r](x, t) (1) 



with 



F[ 



i)]=-0(l + 0')^ G[0(x,t)] = (l + 0^) (2) 



where (/> is a field, 7y is a white Gaussian noise with < 
ri{x, t) >= and < ry(a;, t)'q{x' , t') >= a'^S{x-x')S{t-f), 
and the Stratonovich interpretation is considered pOl. 
Although preliminary simulations seemed to indicate 
that both the NIOT and the NIDT exhibit exponents 
compatible with mean-field like values Q, it has been 
recently elucidated that both transitions belong into the 
kinetic Ising model universality class Q. In particular, 
very extensive simulations and finite size analysis have 
been performed p2[ showing that all the measured ex- 
ponents are compatible with their corresponding Ising 
values. The same conclusion can also be reached by field 
theoretical analysis 0,|3[ ; and it is also extensible to the 
model of Ref. M with both additive and multiplicative 
independent noises. In particular, a naive power counting 
argument permits to conclude that all the terms in eq. 
(0) others than the linear, the cubic, and the Laplacian 
term in the deterministic part, and the constant term 
in the noise, are irrelevant in the renormalization group 
sense. In this way, the relevant part of eq. (H) coincides 
with the Hohenberg and Halperin model A [n3 , and con- 
sequently eq. (y|) is in the kinetic Ising model universality 
class. The irrelevant terms, i.e., the (p^ term in F{<j>) and 
the 4>'^ term in G{4>), play a key role in determining the 
structure of the phase diagram, but do not affect the 
critical behavior. 

In this paper, we face the two previously arisen issues, 
namely, we propose a new candidate with only multi- 
plicative noise exhibiting these two noise transitions, and 
study their critical behavior under a new perspective. We 



will conclude that the phenomena is purely multiplica- 
tive, and therefore the presence of additive noise (i.e. the 
constant term in G{(j))) is not required to generate nei- 
ther the NIOT nor the NIDT. As a byproduct we show 
that the conclusion reached in the previous works |12] 
that both transitions are in the kinetic Ising model uni- 
versality class is not general. We will show that in the 
minimal model we propose, that reproduces the phenom- 
ena under discussion, the transitions belong into a new 
recently proposed universality class: the so called m,ulti- 
plicative noise (MN) universality class |13,nq]. 



II. THE MODEL 

A. Motivation 

A physically intuitive justification of the NIOT has 
been recently proposed [O. It is based on the following 
observation: the short time evolution of the average value 
of the field, (^, of a generic Langevin equation of the form 
(nl) with generic forms of the functionals F and G, is 
easily found to be given by: 

^^A^ = F{< ^>) + ^G(< >)G'(< >) (3) 

when D = {]. 

It is shown in P,0,^.n2| that for the models exhibiting 
a NIOT the previous equation has a linear instability for 
the homogeneous state < (p >= 0; this means that the 
average value of </> grows at initial times, while it would 
decrease monotonously towards in absence of the noise 
term G(< >)G'(< (f) >), i.e., in absence of the multi- 
plicative part of the noise, which is the only one giving 
a non- vanishing contribution to G'{(j))- When the spatial 
coupling is turned on , (i.e., D ^ 0)), it favors neighbor 
sites to take similar values of the field, and, consequently, 
the < (j) >^ solution, induced by the multiplicative 
noise, is stabilized. This is the physical mechanism at 
the origin of the noise-induced phase transition. At this 
point we want to stress the fact that is the multiplicative 
part of the noise the one responsible of the generation of 
the described short-time instability, and therefore of the 
NIOT. 

On the other hand, the NIDT has a more conventional 
interpretation, analogous to standard equilibrium tran- 
sitions: noise (temperature) destroys long range order, 
thus restoring the symmetry of the homogeneous state. 

Inspired by the previous considerations we propose a 
simpler model to analyze the noise-induced transitions 
in which the noise is purely multiplicative. In all the 
previously studied models an additive noise term was also 
present: it is not essential for the phenomena but it does 
change the universality class of the two transitions. 



B. Definition 



Our model is defined by a Langevin equation (M with 
the functionals F and G given by: 

F[c^] = -acj, - cp\ G[0] = 0(1 - a02)i/2 (4) 

where the Gaussian noise, 77, is specified by < rj{x,t) >— 
and 



< ri{x, t)r]{x',t') >= a'^S{x - x')S{t - t'). 



(5) 



and is a positive-definite field, and the Stratonovich 
interpretation is considered. 

In this way, we have a very simple deterministic term 
and a purely multiplicative noise, i.e., with no additive 
part. By setting the parameter a to zero, we get the 
NIOT, therefore the more interesting phenomena un- 
der study is present in the simple multiplicative noise 
Langevin equation studied in Pq ], (that is nothing but 
our model with a = 0). Nevertheless, that model does 
not exhibit a NIDT due to the fact that the order pa- 
rameter keeps on increasing monotonously with increas- 
ing noise amplitude. In order to recover the NIDT a 
number of possible mechanisms can be invoked. Among 
them we have chosen to introduce the term in eq. (0) 
parametrized by a. This "trick" keeps the order param- 
eter from growing indefinitely for large values of a. Let 
us point out that eq. (0) cannot overpass the upper limit 
(f) = cT^I'^ . Observe that while the multiplicative noise is 
essential to generate the NIOT, the NIDT is reproduced 
in a somehow artificial; the NIDT is more naturally in- 
duced by standard additive noise, but it also changes the 
critical behavior. 



C. Renormalization Group results 

Even though we do not present here the details of the 
calculations, it is a straightforward task to conclude from 
a naive power counting analysis that this model belongs 
into the so called multiplicative noise (MN) universality 
class (defined by our model with a = 0, see ||l^). For 
that, it has only to be noticed that a is an irrelevant 
parameter in the renormalization group sense. This con- 
clusion is not altered when introducing diagrammatics 
corrections. Therefore the critical exponents in all the 
possible transitions exhibited by this model are expected 
to be in the MN universality class [l^]. 



D. Mean field results 

We have first performed a mean field analysis of the 
model. For that we consider a given point x and assume 
its nearest neighbors to take a constant undetermined 
value, TO, the expectation value of 4> is calculated as a 
function, / of m, < >= f{m), and the equation is 



closed self-consistently, i.e, < ({> >= m. In figure 1 we 
observe that, for small values oi a, y = /{m) does not 
intersect with the straight line y = m, while they do 
intersect each other for values of a larger than a certain 
critical value, ad- The intersection point is very small 
nearby Cci, grows while increasing the noise amplitude 
up to a point in which starts decaying towards zero for 
large a values (in a extremely slow fashion; see figure 
1). We have shown analytically that in the infinite noise 
limit the only solution is m = for a j^ 0. Therefore, 
in this approximation the model displays a NIOT and a 
NIDT, and both of them are continuous. 

The results for the a = are also displayed in figure 
1(a). Observe that as it was said previously the order 
parameter increases monotonously with a and therefore 
the model exhibits a NIOT but not the NIDT in mean 
field approximation. 
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FIG. 1. (a) Mean field value of < > for different values 
of a, for a = and a — 1. (b) Curves f{in) plotted versus m 
for different values of a and a — 1. The mean field solution 
plotted in (a) is given by the intersection point of the different 
curves with the straight line, y = m (see text). 



E. Numerical results 

We have studied the previous model in both one and 
two dimensions by integrating numerically the associated 
Langevin equation. For that we have used a Runge-Kutta 
type of method (known as the Heun method [|l3,|8| ) . The 
time step is chosen to be St = 10~^, and other parameters 
fixed in the simulations are : a = I, D = 0.2 for the one- 
dimensional simulations, and D = 1 in two dimensions. 
Periodic boundary conditions are employed, and lattice 
sizes up to 1/ = 1000 in d = 1, and L = 40 in d = 2. The 
main results we find are as follows: 

1- We find numerically the two noise-induced transi- 
tions in both d = 1 and d = 2 for a — 1 (for a = we 
find only the NIOT; in the ordered phase the order pa- 
rameter grows monotonously with increasing a and there- 
fore there is no NIDT). In figure 2 the order parameter, 
< (j) > is plotted as a function of the noise amplitude. 



a, for d = 1 and d = 2 (with a = 1). Observe that due 
to the high values of a the finite size effects are much 
more pronounced in the NIDT than in the NIOT, and 
therefore much larger systems and better statistics are 
required to determine accurately the critical point and 
the associated critical behavior of the NIDT. 

Let us stress the fact that we find a pair of phase tran- 
sitions even in d = 1. This is not expected to be the case 
in any of the previously studied models due to the fact 
that in all of them the constant additive part of the noise 
"cohabits" with the purely multiplicative term, and that 
changes dramatically the critical behavior, rendering the 
system Ising-like, and as is well established, systems of 
that nature cannot get ordered in d = 1. 
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FIG. 2. Stationary expectation value of as a function of 
the noise amplitude a in one and two dimensions. 

2- We have performed a detailed numerical analysis of 
the NIOT in one and two dimensions. We have computed 
the critical exponents defined by: 
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In order to locate accurately the critical point and de- 
termine the exponent /x in d = 1 we looked for the value 
of (Jci for which the best linear fit is got in a plot of 
log{< 4> >) versus log[a — ctci)- The exponent ii is given 
by the slope of such a plot with the optimized value of 
CTci- We find: CTcI = 1-80 ± 0.02 with /u = 0.8 ± 0.1 in 
d= 1. 

Afterwards, we determine 9 from the asymptotic decay 
of < 0(i) > right at the critical point, and get 9 = 1.0 ± 
0.1 in d = 1. Next we fix the noise amplitude to its 
critical value and perform simulations changing the value 
of the linear term a. By doing that we find: /3i = 1.50 ± 
0.05 and /32 = 1.7 ±0.2. 

Let us point out that the exponents 9, (3i and (32 are 
in very good agreement with their corresponding values 
in the one-dimensional multiplicative noise universality 
class |lq,f9[ , and the predicted scaling relations hold. In 
particular, we are tempted to conclude from our numer- 
ical simulations that /3i = /32 = 3/2 and 9 = 1 in d = 1 
and, in this way, all the exponents in the MN universality 



class would take rational values. On the other hand, /i is 
a new exponent that had not been previously determined 
for the MN universality class pO| |. 

Ind = 2 concluding results are much tougher to obtain. 

We have verified that the behavior of the computed 
magnitudes at the critical point is compatible with expo- 
nents n = 1.03 ± 0.05, e = 1.75 ± 0.30, P = 1.14 ± 0.05, 
and P2 = 1-70 ± 0.10 |ll,|ll). A more detailed discus- 
sion of this point will be presented elsewhere pit] . Up to 
our knowledge this is the first computation of the critical 
exponents of MN in d — 2. 

3- We have also performed a study of the NIDT in 
d— I and d = 2. As was said previously, finite-size effects 
are much more severe in this transition than they are in 
the NIOT. For that reason, a finite-size-scaling analysis 
has to be performed in order to elucidate its universal- 
ity class. This analysis will be published elsewhere pi[ . 
In any case, based on both theoretical and preliminary 
numerical analysis we believe that these transitions also 
belong into the MN universality class P2|. In particu- 
lar, the fact that the NIDT is present in d = 1 excludes 
the possibility of Ising-like behavior, and relevance argu- 
ments strongly support the previous hypothesis. 



III. CONCLUSIONS 

We have introduced a very simple stochastic model 
exhibiting two noise-induced transitions. We have shown 
that the noise-induced ordering transition (NIOT) is es- 
sentially originated by a purely multiplicative noise, and 
have found it in d = 1 (where it has not been previ- 
ously reported) and d — 2 (where we measure the MN 
exponents for the first time). That shows that the mul- 
tiplicative noise is capable to originate either a NIOT 
or a NIDT even in that low dimension in which phase 
transitions are rare. 

We have verified both analytically and numerically 
that NIOT and NIDT of our model belong into the re- 
cently elucidated multiplicative-noise (MN) universality 
class [^. The result, obtained in previous works, that 
the NIOT and the NIDT belong into the kinetic Ising 
model universality class is not general; that result derives 
from the fact that an additive noise term is also present 
in those models, and that type of term, unnecessary as 
we have shown to generate the phenomena under study, 
does change the system critical behavior, rendering the 
universality class model dependent. 
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We have not succeed so far in relating /i to other known 
exponents using scaling arguments. 
Walter Genovese and M.A. Mufioz. To be published. 
We cannot rigorously exclude the possibility of the NIDT 
being discontinuous (first order). 



